Fluctuation-Dissipation Theorem in Isolated Quantum Dipolar Bosons After a Quench 
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We examine the validity of fluctuation-dissipation relations in isolated quantum systems taken out of equilib- 
rium by a sudden quench. We focus on the dynamics of trapped hard-core bosons in one-dimensional lattices 
with dipolar interactions whose strength is changed during the quench. We find that fluctuation-dissipation re- 
lations hold if the system is nonintegrable after the quench. They also hold if the system is integrable after the 
quench if the initial state is an equilibrium state of a nonintegrable Hamiltonian. However, they fail if the system 
is integrable both before and after quenching. 



pure state \<j>mi) that is n °t an eigenstate of the Hamiltonian 
H (H\\jf a ) = E a \yf a )) (as in Ref. U), then the infinite-time 
average of the evolution of the observable O can be written as 



The fluctuation-dissipation theorem (FDT) OHM] is a funda- 
mental relation in statistical mechanics which states that typ- 
ical deviations from the equilibrium state caused by an exter- 
nal perturbation (within the linear response regime) dissipate 
in time in the same way as random fluctuations. The theorem 
applies to both classical and quantum systems as long as they 
are in thermal equilibrium. In general, for out-of-equilibrium 
systems, fluctuation-dissipation relations are not satisfied. In 
particular, if a system is isolated, it is not clear whether once 
taken far from equilibrium fluctuation-dissipation relations 
will apply at any later time. Studies of integrable models 
such as a Luttinger liquid [4] and the transverse Ising chain 
15[] have shown that the use of fluctuation-dissipation relations 
to define temperature leads to values of the temperature that 
depend on the momentum mode and/or the frequency being 
considered. More recently, Essler et al. 16fl have shown that 
for a subsystem of an isolated infinite system, the basic form 
of the FDT holds, and that the same ensemble that describes 
the static properties also describes the dynamics. 

The question of the applicability of the FDT is of relevance 
to cold atomic systems. In particular, whether it applies to 
few-body observables, such as the momentum distribution 
function, when measured for the entire system. Ultracold 
quantum gases have attracted much attention recently as a 
means for studying equilibrium and nonequilibrium properties 
of model Hamiltonians of interest, and realizing novel quan- 
tum phases 17l|8[]. In such systems, the constituents (bosons 
and/or fermions) and their interactions, as well as the effec- 
tive dimensionality of the region where they are trapped, are 
highly controllable. They are also well isolated from the envi- 
ronment by the ultra-high vacuum in which they are confined 
by means of conservative (magnetic and/or optical) fields. Un- 
der those conditions, one can consider the evolution of the 
system after an applied perturbation to be, to a good approxi- 
mation, unitary |9fl . 

In that context, the description of observables after relax- 
ation (whenever relaxation to a time-independent value oc- 
curs) has been intensively explored in the recent literature 
111 Oil . This is because, for isolated quantum systems out of 
equilibrium, it is not apparent that thermalization can take 
place. For example, if the system is prepared in an initial 
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where c a = (y/«| </>;„;), O aa = (Va|0|Va), and we have as- 
sumed that the spectrum is nondegenerate. The outcome of the 
infinite-time average can be thought of as the prediction of an 
ensemble whose density matrix is diagonal in the eigenenergy 
basis (|c a | 2 being the diagonal matrix elements) ill ill . Equa- 
tion (Q} shows that <9diag depends on the initial state through 
the c a 's (there is an exponentially large number of them) 
while the thermal predictions depend only on the total energy, 
(0i n j|//|0i n j), i.e., they need not agree. 

The lack of thermalization of some observables, in the spe- 
cific case of quasi-one-dimensional geometries close to an 
integrable point, was seen in experiments fll2|l . and, at in- 
tegrability, confirmed in computational 013l4l5ll and analyti- 
cal jlrjll calculations. Away from integrability, computational 
studies have shown that few-body observables thermalize in 
general III U ll7M19ll . which can be understood in terms of 
the eigenstate thermalization hypothesis (ETH) 1 1 lll20l \2 ill . 
We note that the nonintegrable systems studied computation- 
ally belong to two main classes of lattice models: (i) spin- 
polarized fermions, hard-core bosons, and spin models with 
short ra nge (ne arest and next nearest neighbor) interactions 
I II ESS, El, and (ii) the Bose-Hubbard model Cl- 
in this Letter, we go beyond these studies and report re- 
sults that indicate that fluctuation-dissipation relations are also 
valid in generic isolated quantum systems after relaxation, 
while they fail at integrability. For that, we use exact diag- 
onalization and study a third class of lattice models, hard- 
core bosons with dipolar interactions in one-dimension JZj. 
The latter are of special interest as they describe experiments 
with quantum gases of magnetic atoms trapped in optical 
lattices J24ll as well as ground state polar molecules 12511 . 
Rydberg-excited alkali atoms [26Q and laser-cooled ions 027ll 
may soon provide alternative realizations of correlated sys- 
tems with dipolar interactions. The effect of having interac- 
tions beyond next-nearest neighbors in the dynamics and de- 



scription of isolated quantum systems after relaxation is an 
important and open question that we address here. 
The model Hamiltonian for those systems can be written as 
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where b ■ (b •) creates (annihilates) a hard-core boson (£> • = 

b; = 0) at site j, and « = bib: is the number operator. / is the 
hopping amplitude, V the strength of the dipolar interaction, 
g the strength of the confining potential, Xj the distance of 
site j from the center of the trap, and L the number of lattice 
sites (the total number of bosons, p, is always chosen to be 
p = L/3), We set J = 1 (unit of energy throughout this paper), 
fi = ks = 1, use open boundary conditions, and work in the 
subspace with even parity under reflection. 

We focus on testing a fluctuation-dissipation relation after a 
quench for experimentally relevant observables, namely, site 
and momentum occupations (results for the density-density 
structure factor are presented in Ref. |]28|). A scenario under 
which FDT holds in isolated systems out of equilibrium was 
put forward by one of us in Ref. IJ29I1 . There, it was shown 
that after a quantum or thermal fluctuation (assumed to oc- 
cur at time t' , which was treated as a uniformly distributed 
random variable), it is overwhelmingly likely that 0,i± t = 
Cmuc(t)O t >, where O, = (6{t)) H. Formally, C Huc (f) is re- 
lated to the second moments of a probability distribution for 
Ot, Cfiuc(0 = O t+t nO t ii /(O t n) 2 , where infinite-time averages 
have been taken with respect to t". Using Eq. (fl3, and assum- 
ing that no degeneracies occur in the many-body spectrum or 
that they are unimportant, C$i uc (t) can be written as 
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where the proportionality constant is such that Cpi uc (0) = 
1 13 ill - The correlation function in Eq. (01 explicitly depends 
on the initial state through c a . 

Assuming that eigenstate thermalization occurs in the 
Hamiltonian of interest, the matrix elements of O in the en- 
ergy eigenstate basis can be written as: 



O ap = £l(E)S ap +e- s ^ 2 f(E,co)R ap , 



(4) 



where E = A 



(E a +Ep), 0) =E a —Ep, S(E) is the thermo- 
dynamic entropy at energy E, e s ^> = £X a &(E —E a ), £l(E) 
and f(E,co) are smooth functions of their arguments, andR a n 
is a random variable (e.g., with zero mean and unit variance). 
This is consistent with quantum chaos theory and is presum- 
ably valid for a wide range of circumstances 11291 13211 . From 
Eq. ©, it follows straightforwardly that 



Cf1uc(0 ~ CAppr(O) 

where we have defined 
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where again the proportionality constant is such that 
CAppr(O) = 1 0311 . We see from Eq. (0 that Cfi uc (0 does not 
depend on the details of the initial state, in the same way that 
observables in the diagonal ensemble do not depend on such 
details. 

We can then compare this result to how a typical deviation 
from thermal equilibrium (used to describe observables in the 
nonequilibrium system after relaxation) caused by an external 
perturbation "dissipates" in time. Assuming that the perturba- 
tion is small (linear response regime), and that it is applied at 
time t = 0, Coiss(0' defined via O t = Cbiss{jQ0Theimai> can be 
calculated through Kubo's formula as J29l 13411 
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where again, the proportionality constant is such that 
CDiss(O) = 1. Using Eq. ©, one finds that 
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where the last similarity is valid if the width of f(E, Co) 112 811 
is of the order of, or smaller than, the temperature. The re- 
sults in Eqs. © and ([H) suggest that FDT may hold in isolated 
quantum systems out-of-equilibrium under very general con- 
ditions. 

In what follows, we study dipolar systems out of equi- 
librium and test whether their dynamics is consistent with 
the scenario above. This is a first step towards understand- 
ing the relevance of FDT, and of the specific scenario pro- 
posed in Ref. Q29D, to experiments with nonequilibrium ul- 
tracold quantum gases. The dynamics are studied after sud- 
den quenches, for which the initial pure state |0j n ;} is se- 
lected to be an eigenstate of (O for V = Vj n j and g = g; n ; 
(Hj n j), and the evolution is studied under Hfi n (V = Va n and 



g fin ), i.e., \<j)(t)) 



e &at \fyna). We consider the follow- 
ing three types of quenches: type (i) {Vi n ; = 0, gmi = y}^ 
{Vfin = 0, gfi n = y/10} (integrable to integrable), type (ii) 
{Vini = 8, gini = y}^ {Vfin = 0, g fi n = j] (nonintegrable to 
integrable), and type (iii) {V ini = 8, g ini = y}^ {V fin = 2, 
gfin = 7} (nonintegrable to nonintegrable). We choose y such 
that yxj = jx\ = 4, which ensures a (nearly) vanishing density 
at the edges of the lattice in the ground state. The initial state 
for different quenches, which need not be the ground state of 
Him, is selected such thatiitot = (0im|^fin|0ini) corresponds to 
the energy of a canonical ensemble with temperature T = 5, 
i.e., such that^ = Tr{e- fI ^ T Hf la }/Tr{e- a ^ T }. 

In Fig.Q] we show results for Cp\ uc (t ), Coiss( f ), an d Ca PP i-(/ ) 
when the observable of interest is the occupation of the site 
in the center of the system n ;= L/2 (qualitatively similar results 
were obtained for other site occupations, for momenta occu- 
pations, and for the density-density structure factor lt28IO . The 
results are obtained for the three different quench types men- 
tioned above and are shown for L — 15 and 18. For quench 
type (i), we find that none of the three correlation functions 
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FIG. 2. Normalized variance of CpiucM -QmssW an d Cpi uc (f) — 
CAp pr (0 vs the system size for the three quenches explained in the 
text [identified by Q (i), Q (ii), and Q (hi)]. The observable is the 
occupation of the site at the center of the system. The variances are 
calculated for 2000 points between t = and 100, and normalized by 
the average variance of the two correlation functions for which the 
differences are computed. 



FIG. 1. Correlation functions, CpiucM, CpjissWi an dCA PP r(0 wnen 
the observable is W;=l/2 vs ti me *■ Results are shown for the three 
quenches (i)-(iii) (from top to bottom, respectively) explained in the 
text, and for L = 15 (left panels) and 18 (right panels). Results for 
L = 12 are presented in Ref. 112811 . The insets show normalized his- 
tograms of Cpi uc (f) (red filled bars) and Cd; ss (/) (blue empty bars) 
calculated for 2000 data points between / = and 100. 



agree with each other and that the agreement does not im- 
prove with increasing system size [see Figs. Q2a) and[TJb)]. 
There are also large time fluctuations, characteristic of the in- 
tegrable nature of the final Hamiltonian II35I1 . We quantify 
these fluctuations by plotting the histograms of Cfi uc (0 and 
Crass (0 f or an extended period of time in the insets. We find 
the histograms to be broad functions for quenches (i) and (ii) 
[see Figs, fflftffld)]. 

Remarkably, in quenches type (ii) [Figs. [TJc) and |TJd)], 
which also have a final Hamiltonian that is integrable, Cpi uc (f ) 
and Crass (0 are very similar to each other and their differ- 
ences decrease with increasing the system size. This indicates 
that the FDT holds. At the same time, we find differences 
between fluctuation/dissipation correlations and Ca PP i-(0> m ~ 
dicating that the agreement between Cfi uc if ) and Crass ( f ) does 
not imply that Eq. (0|i is valid. These observations can be un- 
derstood if the initial state provides an unbiased sampling of 
the eigenstates of the final Hamiltonian. In that case, even 
though eigenstate thermalization does not occur, thermaliza- 
tion can take place 13611 and this results in the applicability 
of FDT. In quenches type (ii), such an unbiased sampling oc- 
curs because of the nonintegrability of the initial Hamiltonian, 
whose eigenstates are random superpositions of eigenstates of 
the final integrable Hamiltonian with close energies 113611 . 

For quenches type (iii) [Figs. [T]e) and [Iff)], on the other 
hand, we find that not only Cpiucif) an< ^ Crass (0 are very close 



to each other, but also C^-pmit) is very close to both of them, 
and that the differences between the three decrease with in- 
creasing the system size. Therefore, our results are consistent 
with the system exhibiting eigenstate thermalization, which 
means that the assumptions made in Eq. (0]i are valid, and the 
applicability of the FDT follows. We have studied the ther- 



malization properties of this system in more details 112 811 and 
found that they are qualitatively similar to those in models 
with short-range interactions IU7II . Furthermore, for type (iii) 
quenches, one can see that time fluctuations are strongly sup- 
pressed when compared to those in quenches type (i) and type 
(ii) [better seen in the insets of Fig.[TIe)-(f)], which is a result 
of the nonintegrable nature of the final Hamiltonian 0291 137 ] . 
To quantify the differences between the three correlation 
functions and explore their dependence on the system size 
for each quench type, we calculate the normalized variances 
of CfiucO) - Crass (*) and Cfi uc (0 -C A p P r(0> where the nor- 
malization factor is the average variance of the two func- 
tions for which the differences are calculated, e.g., Var(CFiuc — 
CDissV^rVarfCFiuc) -f-Var(Cra ss )]. In Fig. [2j we show these 
quantities for the three quench types vs L. For quench type 
(i), the variances exhibit a tendency to saturate to a non-zero 
value as L increases, which indicates that Cp\ uc (t) and Crass(0> 
as well as Cp\ uc (t) and Ca PP i-(0' ma y remain different in the 
thermodynamics limit. This is consistent with the findings 
in Refs. J4i |5|] where it has been shown that conventional 
fluctuation-dissipation relations with a unique temperature do 
not hold in integrable systems. For quench type (ii), we see 
that the variance of Cfi uc (t) — Coiss (t) decreases with increas- 
ing system size and becomes very small already for L = 1 8, in- 
dicating that Cpi uc (t) and Crass (*) agree in the thermodynamic 
limit. The variance of Cp[ uc (t) — Ca PP i(0> on me other hand, 
exhibits a more erratic behavior, and it is not apparent whether 
it will vanish in the thermodynamic limit. For quench type 
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FIG. 3. Absolute value of the off-diagonal matrix elements of n j = lI2 
and «£=o m me eigenenergy basis, in a narrow energy window 
around E = Et \ (with a width of 0.1), vs the eigenenergy differ- 
ence, (0 = E a — En. Results are shown for L = 15 (left panels) and 
L = 18 (right panels). (a),(b) correspond to the final Hamiltonian 
in quenches type (ii), and (c),(d) correspond to the final Hamilto- 
nian in quenches type (iii). The green (light gray) symbols are the 
off-diagonal matrix elements of i\ = q, and the black ones of fij= L /2- 
Because of the overwhelmingly less number of data points for n ; = l/2 
and their sparsity in the ranges shown for the integrable case, in 
(a),(b), we have increased the size of the symbols for rc.- =z / 2 by a 
factor of 20 relative to those for «j = o, making them distinguishable 
from each other. To increase the resolution of the distribution of val- 
ues in the case of L = 18, where a very large number of data points 
exists, we plot only 1 out of every 10 points for «£ = o in (b), and for 
both observables in (d). Lines are running averages for n,t = o with 
a subset length of 50 for L = 15 and 200 for L = 18. Insets show 
the histograms of the relative differences between the «t = o data and 
running averages with subset sizes of 1000 for L = 15 and 10000 for 
L = 18. The relative difference is defined as (|O a R| — /avg)//avg, 
where / avg is the running average. 



(iii), the relative differences between Cfi U c(0' Cjy{ ss (t), and 
CAppr(f ) exhibit a fast decline with increasing the system size, 
indicating that all three will agree in the thermodynamic limit. 
These results strongly suggest that the FDT is applicable in the 
thermodynamic limit for quenches in which the final system 
is nonintegrable, as well as after quenches from nonintegrable 
to integrable systems, even though the ETH does not hold in 
the latter. 

In order to gain an understanding of why FDT fails or ap- 
plies depending on the nature of the final Hamiltonian, we 
explore to which extent Eq. (HJi describes the behavior of the 
matrix elements of few -body observables in the nonintegrable 
case, and in which way it breaks down at integrability. In 
Fig. [3] we plot the off-diagonal elements of two observables, 
iij =L i 2 , and the zero-momentum occupation number, n^. = o (see 



Ref. B28I1 for further results) vs the eigenenergy differences 
(CO) in a narrow energy window around E = E m . Results 
are shown for matrix elements in the eigenstates of the final 
Hamiltonians in quenches type (ii) and type (iii) 13811 . The 
off-diagonal matrix elements of both observables in the eigen- 
states of the integrable Hamiltonian [Fig. [3ja)-[3jb)] can be 
seen to exhibit a qualitatively different behavior from those 
in the nonintegrable one. In the integrable Hamiltonian, the 
matrix elements exhibit extremely large fluctuations. In addi- 
tion, a very large fraction of those elements (larger for Hj=l/2 
than for H£ = o) have vanishing values. This makes any def- 
inition of a smooth function, f(E,co), meaningless. Those 
results contrast the ones obtained in the nonintegrable case, 
where the fluctuations of the matrix elements have a different 
nature, and we do not find a large fraction of vanishing ones. 
To see that more clearly for ii^q (the better behaved of the 
two observables), in the insets of Fig.|3]we show the normal- 
ized histograms of the relative differences between the matrix 
elements for tik=o and a "smooth" function, defined as the run- 
ning average of those elements over a large enough group of 
them (examples of the running averages are presented in the 
main panels). For the integrable system, we find that the his- 
tograms are not compatible with the uniform distribution pos- 
tulated in Eq. @, as a very sharp peak develops at — 1 for both 
system sizes. That peak becomes increasingly sharper with in- 
creasing system size, reflecting the increasingly large number 
of vanishing off-diagonal matrix elements in those systems. 
For the nonintegrable Hamiltonian, on the other hand, the his- 
tograms are closer to a uniform distribution, though they do 
have a small finite slope. 

In summary, studying the dynamics of an experimentally 
relevant model of trapped hard-core bosons with dipolar in- 
teractions, we have found indications that the fluctuation- 
dissipation theorem is applicable to the properties of few-body 
observables in nonintegrable isolated quantum systems out of 
equilibrium, and that this follows from the ETH. Furthermore, 
we find indications that the fluctuation-dissipation theorem 
may also apply to integrable systems, for which the ETH is 
not valid, provided that the initial state before the quench is 
an equilibrium state (eigenstate) of a nonintegrable system. 
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EXPERIMENTAL RELEVANCE OF OUR MODEL 

Hamiltonian (2) in the main text provides a microscopic 
description for the dynamics of a gas of bosonic ground 
state molecules such as, e.g., LiCs molecules (dipole moment 
flWx ~ 5.6 Debye), confined transversely (longitudinally) by 
a two-dimensional (one-dimensional) optical lattice with fre- 
quency a>j_ (0)||), with 0)x 3> ©|| ■ The molecules are polar- 
ized in the transverse direction by an external electric field of 
strength F and are confined to the lowest band of the ID lat- 
tice, provided o>i > max{V,J,T}. Here, V = d 2 /{AkeqoI), 
with d < c/ m ax the dipole moment induced by F, an the lat- 
tice spacing in ID and £o the vacuum permittivity. The hard- 
core condition is obtained by requiring that molecules are 
trapped with a low-density n, such that the initial system 
has no doubly occupied sites [39]. The additional condition 
re -1 ' 2 ^$> (djLa/ftWll) — 360nm ensures collisional stabil- 
ity 13911 . The model in Eq. (2) of the main text can also 
describe the dynamics of a gas of strongly magnetic atoms 
such as Dy (dipole moment d = lOjits, with jig Bohr's mag- 
neton) or Er (d = 7jUb). In this case, the hard-core constraint 
is achieved by means of magnetic tuning of the short-range 
scattering length using Feshbach resonances [40], while J/V 
decreases exponentially with increasing the depth of the ID 
lattice Jl. 

THERMALIZATION 

Despite the presence of interactions that have a power-law 
decay with distance, we find that the behavior of eigenstate ex- 
pectation values of few-body observables, as well as thermal- 
ization properties of the systems described by Hamiltonian (2) 
of the main text, are qualitatively similar to those already seen 
in models with short range (nearest and next nearest neighbor) 
interactions 11171 kin . 

In order to show that this is indeed the case, here we study 
the difference between the results of the diagonal ensemble for 
the few-body observables studied in main text, namely, tij and 
Ufa as well as for the density-density structure factor, iVj (not 
studied in the main text), and the results of the microcanon- 
ical ensemble. The momentum distribution function and the 



density-density structure factor are defined as 

W&'-^tfc, N k =\l,eW- m \n m . (9) 



"* 



' l,m 



They are the Fourier transforms of the one-particle and 
density-density correlation matrices, respectively. Since we 
work at fixed number of particles, (A^=o) = P 2 /L, so we set it 
to zero without any loss of generality. These observables can 
be studied in ultracold gases experiments. 

We define the microcanonical average for an observable O 
as <9 m icro = -jr- La O aa . Here, jV^e is the number of states in 
the microcanonical window, which is centered around E tot and 
has a width of AE l42ll . We average the results over several 



close values of AE for each system size to ensure that they are 
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FIG. 4. (a)-(c) Normalized differences between the diagonal and mi- 
crocanonical ensemble averages for the three observables, rij, n^, and 
iVj [see Eq. UOH vs L. Different plots in each panel correspond to 
the three quench types [identified by Q (i), Q (ii), and Q (iii)]. (d)-(f) 
Maximum of the normalized differences between the microcanonical 
average and each of the eigenstate expectation values in the micro- 
canonical energy window [see Eq. i ll IH vs L. 




0.6 

s 

II 

e~" 0.3 




R *0.5 



0.4 

K 

II 

fe; 0.2 




FIG. 5. Eigenstate expectation values values of n;^ L n, %=o> an ^ 
ivjt = ji over the entire spectrum vs the scaled eigenenergies for the 
noninteracting (a)-(c), and interacting system (d)-(f). We show re- 
sults for L = 15 as the green (light colored) points, and for L = 18 
as the black (dark) ones in each panel. The vertical lines show the 
location of E tot /L for L = 18, which is also roughly equal to that of 
L=15. 



robust against small changes in AE. The values for AE are in 
[0.2 - 0.25] for L = 12, and [0.1 -0.15] fori = 15 and 18. 

We then compute the normalized difference between the 
predictions of the diagonal and the microcanonical ensembles, 






0.5 
E /L 



AO 



EflQdiagffl-Qmtaoffll 



(10) 



where £ stands for either the site index or the momentum index 
depending of the observable. 

In Fig. |4|a)-(c), we show results for AO vs L for the three 
types of quenches described in the main text. We find that for 
quench type (i) (between two integrable systems), AO is gen- 
erally larger in comparison to AO for the other two quenches, 
and that this feature remains apparent as the system size in- 
creases to L — 18. Previous studies, involving the same inte- 
grable model without Il3l4l5ll and with 11411 the trapping po- 
tential and much larger system sizes, have found strong indi- 
cations that AO (where O was either the density or the mo- 
mentum distribution function) remains finite in the thermody- 
namic limit so that the system does not thermalize. On the 
other hand, trends for quench type (iii) suggest that the differ- 
ences vanish by increasing the system size, i.e., that the non- 
integrable dipolar system thermalizes in the thermodynamic 
limit. The trend is less clear for quench type (ii) as we see 
that Anj monotonically decreases by increasing L, whereas 
An/t and AiVj appear to saturate by increasing L from L = 15 
to L = 18. Clearly, larger system sizes are required to con- 
firm whether thermalization takes place for such an integrable 
system after a quench from a nonintegrable one as argued in 
Ref. Q. 



As mentioned in the main text, for a generic (nonintegrable) 
system, thermalization can be understood through the ETH. 
Here, we examine the validity of the ETH for each case con- 
sidered by calculating the normalized differences between the 
observable in each eigenstate and the microcanonical average, 



AO c 



L/|Oaa(Q-Omicro(/)| 
L/<2micro(0 



(ID 



and taking the maximal difference within the microcanoni- 
cal window, [A(9 aa ] max = Max{AO a a]AE- This is a measure 
of how widely the eigenstate expectation values are spread 
in the microcanonical window. In Fig. Etd)-|4tf), we show 
this quantity for our three observables vs L. They are seen 
to consistently decrease with increasing system size for the fi- 
nal Hamiltonian in quench type (iii) while they are seem to 
saturate to relatively large values for the final Hamiltonians 
in quenches type (i) and type (ii). This is an indication that 
eigenstate thermalization occurs in the former case while it 
fails in the latter ones. 

The validity, or failure, of ETH can be perhaps more easily 
seen by plotting the eigenstate expectation value of observ- 
ables vs the eigenenergies, as shown in Fig. For the inte- 
grable system with V = [Fig.[3a)-|3c)], eigenstate expecta- 
tion values exhibit large fluctuations inside the microcanoni- 
cal energy window for both L = 15 and 18, and the width of 
the region where the values are scattered around E tot does not 
decrease with increasing the system size, i.e., eigenstate ther- 
malization does not occur. This is different from what hap- 
pens in the nonintegrable case [Fig. [5jd)-|5tf)], where each of 
the eigenstate expectation values of an observable inside a nar- 
row energy window around E tot approaches the microcanoni- 
cal average as the system size is increased. This is apparent as 
the width of the region where the values reside around E tot de- 
creases with increasing system size, and presumably vanishes 
in the thermodynamic limit. 

ENERGY SCALES 

The short-time evolution of the correlation functions is set 
by the width of f(E,a>) as a function of (0, which we de- 
note as W. Note that, as discussed in the main text, / is not 
well defined for integrable systems and, even for the nonin- 
tegrable case for which Eq. (4) of the main text is seen to 
better describe the data as the system size increases, one can- 
not disentangle / from the random function R a p with merely 
the knowledge of the off-diagonal values. Therefore, we esti- 
mate W for another closely related function, / cg (a>), which is 
obtained by coarse-graining the off-diagonal values. We then 
calculate the width using 



W 



2f~\f cg (co)\ 2 dco 

l/cg(0)| 2 ' 



(12) 



Examples of / cg (o)) are depicted as lines in the right panels 
of Fig. [7] We choose different bin sizes for different systems. 
They are, 50 x 8 for L = 12, 300 x 8 for L = 15, and 700 x 8 
forL =18, where 8 is the average level spacing. 
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FIG. 6. Various energy scales in the model vs L. Top panels show 
the estimates for the width of the off-diagonal function, f(E, ft)) in 
the ft) space, (d) shows the uncertainty in the total energy, and (e)-(f) 
show the average level spacings (see text). 



In Fig. [6ja)-[6jc), we plot W versus L in the three quenches 
and for the three observables considered. For « j=l/2 and n k=o, 
W varies non-monotonically with L in most cases, while for 
Nit =n , it is seen to monotonically decrease with increasing L 
and possibly saturate to a finite value for larger system sizes. 
Larger system sizes are needed to understand the behavior of 
W in macroscopic systems. Regardless, the values obtained 
here can be used in each case to estimate the time scale for 
the initial decrease of the fluctuation or dissipation correlation 
functions. 

In Fig. |6jd), we also show the quantum uncertainty of the 
energy, A. We find that they are ^(lO ) for the systems stud- 
ied here, and increase sl owly with system size, as expected 
from the analysis in Ref. 1 1111 . 

The time scale for recurrences in the correlation functions is 



set by the average level spacing, 8 (estimated by e s ^ E 'E Il29l0 . 
which, as expected, is found to decrease exponentially with in- 
creasing the system size from ^(10 _1 ) forL = 12 to ^(10~ 3 ) 
forL= 18 [see Fig. [(Je)]. Therefore, such a time scale for typ- 
ical system sizes explored in experiments would be much too 
large to have any relevance. This is also true if one calculates 
<5eff (= v t- i-t ) II22D' snown m Fig-IUf), which represents the 



Y,a \ c a\ 



effective level spacing between the eigenstates participating in 
the diagonal ensemble. 

OTHER OBSERVABLES AND/OR SYSTEM SIZE 



In the main text, we show the fluctuation and dissipation 
correlation functions of « j= l/2 f° r me two largest system 
sizes, L = 15 and 18. For completeness, in Fig. [7] we show 
the same quantities, as well as the corresponding off-diagonal 
elements of the two observables shown in Fig. 3 of the main 
text, for the smallest system we have studied, L — 12. Because 
of the smaller size of the Hilbert space in comparison to the 
other clusters, some of the trends seen in Fig. 1 of the main 
text are not so clear in the case of L = 12. However, the sup- 
pression in the fluctuations of the correlation functions for the 
nonintegrable case [Fig. |7jc)] is apparent in the correspond- 
ing histogram. Other features, such as the dramatic change 
in the behavior of the off-diagonal elements of «/=£/2 when 
the integrability-breaking interaction is introduced, is already 
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FIG. 7. Left panels (Right panels): the same as Fig. 1 (Fig. 3) 
of the main text, except that L = 12, and that the lines in the right 
panels are coarse-grained functions of the off-diagonal values used 
to calculate the width of data in ft) using Eq. (l2\ . In the right panels, 
solid (dashed) lines correspond to Mjt=o («/=l/2)- The bin size for 
coarse-graining is 50 x S for this system size. 

In Figs.[8]and|9] we show results for the fluctuation and dis- 
sipation correlation functions of «^ = o and A^ =7r for the three 
system sizes and quench types considered in the main text. 
The contrast between the results for different quenches can 
be seen to be similar to the one in Fig. 1 of the main text 
for Hj =L i 2 - The main difference between the results for «^ = q 
and Nk=Tt when compared to those for « j= l/2 i s that the for- 
mer two exhibit smaller fluctuations with increasing system 
size than the latter one. This is to be expected as the presence 
of the harmonic trap, which breaks translational symmetry, 
produces a larger number of nonzero off-diagonal matrix ele- 
ments of «^ = o and Nk =Jl in the integrable regime than of nj =L i 2 
(see Fig. 3 of the main text). Still, those fluctuations [see the 
histograms in the insets of Figs.|8la)-[8jf) and Figs.[9ja)-|9jf)] 
can be seen to be much stronger than in the nonintegrable 
case [see the histograms in the insets in Figs. [HJg)-[8ji) and 

Figs.iagHSi)]. 
The results for the scaling of the variances between 

C F iuc(0 -C D iss(0 and C F i uc (0 -C A p P r(0 for n k=0 and N k=n 
are presented in Fig. [10] They are also consistent with the 
conclusions extracted from the scaling of those differences for 

n j=L/2- 

Finally, in Fig. QT1 we show results for the off-diagonal ma- 
trix elements of Nk =Jl within AE of E tot . Those results are the 
equivalent of the ones presented in Fig. 3 of the main text 
for nj =L i2 and «,t = o. That figure shows that the conclusions 
drawn for the latter two observables in the main text are also 
applicable to Nk =K . 




FIG. 8. The same as Fig. 1 of the main text, except that here, the observable is the zero-momentum occupation number, «£=()• We have also 
included the results for L = 12 in the left panels. 
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FIG. 9. The same as Fig. 1 of the main text, except that here, the observable is the k = lt density-density structure factor, Nk=%- We have also 
included the results for L = 12 in the left panels. For quench type (i), we find that the value of N^ =K for L = 12 does not exhibit any dynamics. 
This is reflected in the corresponding fluctuation correlation function in (a), which is zero at all times. For this reason, we have set Cp] uc (t) = 
in that case and are not showing its histogram. This is unique to N/, =n for L = 12 and to the specific parameters chosen for quench type (i). 
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FIG. 11. Same as Fig. 3 of the main text, but for when the observable is Nk =7Z . Here, we show all the data points in all of the cases. Lines are 
running averages with a subset length of 100 for L = 15 and 1000 for L = 18. 



